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Abstract 



This paper studies the controllability problem of a parabolic system of chemotaxis. The local 
exact controllability to trajectories of the system imposed one control force only is obtained by 



I/"") , applying Kakutani's fixed point theorem combined with the null controllability of the associated 

linearized parabolic system. The control function is shown to be in L°°(Q), which is estimated 
, by using the methods of maximal regularity and L p -L q estimates of parabolic equations. 
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a 

1 Introduction and main results 

Let O C M N (N > 1) be a bounded domain with sufficient smooth boundary dQ. Let w be a 
nonempty open subset of $7, and T > 0. We denote Q = f2 X (0, T), E = d£l x (0, T) and 
= wx(0,T). Throughout this paper, we use W s ' q {Q), Wq ,l {Q) &ndC a (Tl) (s, a > 0, 1 < q < oo) 
for the usual Sobolev spaces (e.g., [30]), and set H m {£l) = W m ' 2 (n) for m 6 N. L p (tt) and I?{Q) 
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(1 < p < oo) are the usual Lebesgue function spaces with the norm | • L 
Moreover, let 



and 



IP' 



respectively. 



V\Q) = {y\yeL 2 {^T;H l (a)),d t yeL 2 {^T;H l ^ir)}, 
V 2 (Q) = {y\yeL 2 (0,T;H 2 (n)),d t yeL 2 (Q)}, 

be equipped with their graph norms, where H 1 ^)* denotes the dual space of The duality 

between H 1 ^)* and H l (Q) is denoted by (•, •). 

In this paper, we are concerned with the following controlled parabolic system with state func- 
tions u = u(x, t) and v = v(x, t) : 



d t u = V • (V-u - x«Vv) + lujf in Q, 

dtv = Av — jv + 5u in Q, 

d v u = 0, d v v = on S, 

u(x, 0) = uq(x) v(x,0)=vq(x) x G fi, 



(1.1) 



where dt = d/dt, and d u = d/dv stands for the derivative with respect to the outer normal v of 
dQ, 1^ represents the characteristic function of u, f = f(x,t) is the control function so that \ u f 
is the control force acting from the outside on a portion of the domain Q, uq and i>o are the initial 
values, and x, 7 and 5 are given positive constants. 
A pair of functions (u, v) with 

u G V\Q) n L°°(Q), v G V 2 (Q) n L°°(Q) 

is called a weak solution of (jl.ip if for all if G L 2 (0, T; the following identities hold: 

f (d t u, (p) dt + J J [{Vu - xuVv) ■ Vip + l u fip] dxdt = 0, 
Jo J Jq 

ipdtvdxdt + / / [Vu • V93 + (—yu + 5u)ip] dxdt = 0. 



We write the free system of (jl.ip . that is, in the absence of /, as follows: 



dtu = V • (Vu — xuSJv) 
dtv = Av — jv + Su 
d v u = 0, d v v = 



in Q, 
in Q, 
on S, 



(1.2) 



«(x, 0) = uq(x) f (x, 0) = vq(x) xEfl. 



The system (|1.2p is a prototype chemotaxis system so called Keller-Segel model which describes 
the aggregation process of slime mold resulting from chemotactic attraction. In (jl.2p . u represents 
the density of the cellular slime mold, v is the density of the chemical substance (see [29J). In the 
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last decade, there is a large number of works devoted to the mathematical analysis of the Keller- 
Segel system. Several topics on the Keller-Segel model for chemotaxis such as aggregation, blow-up 
of solutions, and chemotactic collapse, etc., have been concerned and some significant results have 
been achieved from different discipline perspectives. In Horstmann [27] and Hillen and Painter |26j, 
it provides a detailed introduction into the mathematics of the Keller-Segel model for chemotaxis 
with abundant references therein. Here we would mention a few facts about the local and global 
existence of solutions for the Keller-Segel model. Generally speaking, the blow-up of solutions of 
Keller-Segel system in finite or infinite time depends strongly on the space dimension. In 1-d case, 
a finite time blow-up never occur, and the global solution exists and converges to the stationary 
solution as times goes to infinity (see |32|). But the blow-up may occur in finite or infinite time in 
n-dimensional case for n > 3 (see [11[ 128]). For the 2-d case, several thresholds have been found. 
When the mass of the initial data is below some threshold value, the solution exists globally in 
time and its L°°-norm is uniformly bounded for all time. While the mass of the initial data is 
larger than some threshold value, the solution will blow up either in finite or in infinite time (see 
P [201 [35]). 

Due to blow-up feature of solutions of the Keller-Segel model, it is interesting to consider 
some controllability problems. Let (u, v) be a trajectory, i.e., a solution of (j 1.2)) corresponding 
to some initial value (uq,vq). We say that the system (|1.2p is locally exactly controllable to the 
trajectory (u,v) at time T, if there exists a neighborhood O of (uq,vq) such that for any initial 
data (uq,vq) € O, the solution (u,v) of (11. ip driven by some control function / satisfies 

u(x, T) = u(x, T), v(x, T) = v(x, T), for x € a.e., 

where the neighborhood O and the control function space will be specified later. 

In this paper, we suppose that u,v verify the following regularity properties: 

u,v £ L°°(Q), VtJ G L°°(Q) N . (1.3) 

Remark 1.1. The solution (u,v) of the system (|1.2p exists at least locally in time interval [0, T\] 
for Ti < 

T m ax with sufficiently small initial data {uq,vq), where T max is the maximal existence 
time (see |27] and reference therein). If the system is locally exactly controllable, then we can 
drive the state of the system by some control force to a given trajectory at time T <T\ before the 
time T max to avoid blow-up. It is also worth indicating that the reason we consider the local exact 
controllability instead of exact controllability is that the solution may blow up when the mass of 
initial value is larger than some threshold value. 

Remark 1.2. When (u,v) = (0,0), the local exact controllability is reduced to the local null 
controllability. If the system (jl.ip is locally null controllable at time T with some control, then we 
can switch off the control after time T and the system will keep into zero afterwards. 
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This paper is devoted to the local exact controllability of the coupled parabolic system (jl.ip 
via one control. The controllability of parabolic systems of coupled equations attracts intensive 
attention in the last few years. In Barbu [7], it studies the local exact controllability to steady 
states with controls acting on each equation of the system via the same interior domain. This 
could be done by taking it as a direct consequence of the controllability of the scalar parabolic 
equations. It is much more interesting and applicable to consider the controllability of a parabolic 
system with one control force imposed on one equation of the system. Ammar Kdjodia et al. pQ is 
the first work of this kind. They show that the phase-field system is locally exactly controllable to 
the trajectory by one control force. The series of works of Ammar Kdjodia et al. ([21 [HJ [3]), and 
the works of Gonzalez-Burgos et al. ( |22j [10\ [T4"l [T5] ) , have extended such problem to more general 
cases. The survey paper [5] gives a comprehensive introduction to this topic. For more works of 
the controllability of parabolic equations, we also refer to [161 1171 \TE[ [TO] and |12j . 

However, as to our best knowledge, very few results are available to the control problems of the 
system . In Ryu and Yagi [33] , it considers an optimal control problem of the system with 
the control to be distributed on the second equation of (jl.lj) . The present paper can be considered 
as a first work on the controllability of the system (jl.lj) . There are some other kinds of interesting 
control problems for the system (jl.ip . In the system (jl.ip . the chemotactic term — xV • (uVu) 
causes much more mathematical difficulties than the coupled parabolic systems aforementioned. 
The techniques presented in this paper would be useful for other forms of chemotaxis system such 
as the parabolic-elliptic chemotaxis system, and even for other coupled systems like drift-diffusion 
equations from the semiconductor device. 

The idea of obtaining the controllability of (jl.ip is somehow classical: We first establish the null 
controllability of the linearized system and then apply the fixed point theorem. Now we consider 
the null controllability of the linearized system of (jl.lj) . which is written as follows: 

f dty = Ay - V • (By) - V • (oVz) + l u f in Q, 
dtz = Az — 7Z + 5y in Q, 

(1.4) 

d v y = 0, d u z = on S, 

y(x,0) = y (x) z(x, 0) = z (x) x G Q, 

where a G L°°(Q), B G L°°(Q) N with B ■ u = on S, / G L 2 (Q) is the control force, and 
yo,ZQ G L 2 (f}) are given initial data. To study the null controllability of (jl.4p . we are led to 
consider the observability of the adjoint system of (jl .4j) : 

-dt<f> = Acj) + J3V0 + 69 mQ, 
-d t 9 = AO - 7 - V • (aV(j>) in Q, 

(1.5) 

d u (f> = 0,d u 9 = on E, 

<f>(x, T) = T (x),8(x, T) = 6 T (x) ie!l, 
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where 4 >T ,9 T G L 2 (J7). It is well-known that the null controllability of (jl.4p is equivalent to the 
observability inequality for system (jl.5p : 

\<f>(;0)\l + \e(;0)\l<C [[ \<f>\ 2 dxdt 

J J Qui 

for every solution (4>,6) of (11. 5h . However, in order to obtain the input space of L°°(Q), we need 
to establish instead an improved observability inequality of the following 

\<K;0)\l + \0(;0)\l<C ff e 3 2 sa \<P\ 2 dxdt, (1.6) 

J J Qui 

which can be derived from a global Carleman inequality 

jje 2sa (\<l)\ 2 + \e\ 2 ^dxdt<cjj e'l sa \cj)\ 2 dxdt (1.7) 

for every solution ((j>,6) of (|1.5p . Here, in (|1.6p and (jl.7p . C denotes some positive constant 
independent of <fi and 9, a = a(x, t) is a weight function which will be specified precisely in Section 
2, and s is a real number considered as parameter. The basic idea for the inequality (jl.7p comes 
originally from [T3] and [22], where similar inequalities are obtained for some cascaded system and 
parabolic system of phase-field. 

Now we state our first result. 

Theorem 1.1. Let T > 0. For any (yo,£o) £ L 2 (^) x -L 2 (0), there exists a control f G L°°{Q) 
such that the solution (y,z) of system (|1.4|) corresponding to f satisfies (y,z) G V-^Q) x 1 /1 (Q) 
and T) = 0, z(x, T) = for x G almost everywhere. Moreover, the control f satisfies 

\\f\L<e CK (\yo\ 2 + \\zo\\ 2 ) (1.8) 
where C is a positive constant depending only on £1 and cj, and 

k=(1 + Mi, + WBWlJT + i + 1 + ||o||oo + ||B||oo. (1-9) 

The approach used here to obtain the L°°(Q) control is originally from [7] (see also [36]). We 
improve this approach to get the explicit representation of the bound with respect to T by adopting 
some techniques from semigroup theory such as L p -L q estimate and maximal L p -regularity. 

The main result of this paper is the following Theorem 11.21 

Theorem 1.2. Let p > N + 2. Let (u,v) be a trajectory of (|1.2p corresponding to (uo,vo) and 
satisfy (jl.3p . Then, there exists a positive constant c\ independent of T such that for each {uq,vq) 
that satisfies 

l«o - «oL + IN - ^o|| w2(1 _i), P(n) < e- Cl ( 1+T+ T), (no) 
there is a control f 6 L°°(Q) such that system (jl.ip admits a solution (u,v) satisfying 

u g v\Q) n l°°(Q), u g y 2 (Q) n l°°(Q), 

and it(x, T) = u(x, T), u(a;, T) = v(x, T) for x G almost everywhere. 
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We proceed as follows. In next section, Section 2, we give some preliminary results. Section 3 
is devoted to the proof of the Theorem 11.11 The proof of Theorem 11.21 is presented in section 4. 

It is pointed out that throughout the paper, we use C to denote a positive constant that is 
independent of time T in most cases but may be dependent of In the later case we may write 
C(£1,lo) instead of a special specification. 



2 Preliminaries 

In this section, we collect some results that are needed in later sections. These results are particu- 
larly useful in the establishment of the regularity of linear parabolic system and the L°°-estimate 
of controls. 

For p G (1, oo), let A := A p denote the sectorial operator defined by 

A p u := -An, Vm£ D(A p ) := {u G W 2 ' P {Q); d„u\ an = 0} . (2.1) 
Suppose that 7 is a positive constant. 

(i) Let a > and D ((A + 7)°) be the function space endowed with the graph norm. Then 
D ((A + 7)") is a Banach space with the following embedding properties ( |24| p. 39]) 

D({A + ~f) a ) -4 W lj, (Sl) ifa>i (2.2) 

and D((A + j) a ) ^ C 7 (0) if0< 7 <2a--. (2.3) 

P 

(ii) Let {e~ tA } t>0 and je - "^"*" 7 ) } t>Q be the analytic Co-semigroups generated by —A and — (A+ 
7) on L p (r2)(l < p < 00), respectively. By standard Co-semigroup theory, we have f[T3"l [33] ) 



e u\ q < Cm{t)~^ ( v— q >\ u \ p , (2.4) 



and 



(A + j) a e- t{A+ ^ < Ct-- ( p—v>- a \u\ p (2.5) 



N 1 1 In 



_ V n a 1 



for all u € L p (£l),t > and 1 < p < q < 00, where m(t) = min{l, t}. 

(iii) Let a > and 1 < p < 00. Then for any e > there exists a constant C £ depending on $7, e 
and p such that ([281 Lemma 2.1]) 

I {A + 7 )« e -*A v . u \ p < <7 6 i-<H-e | M | p (2.6) 

for all n G D>(Sl),t > 0. 
As a consequence of (|2.4p and (|2.6p . we have 
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(iv) For any e > 0, there exists a constant C £ depending on f2, s and p, such that 

\e~ tA V -u\ < C £ m(tY 
for all u £ L P (Q), t > 0, 1 < p < q < oo. 



-l-e-^fl-il i , 
2 e 2 u 



(2.7) 



(v) (Maximal regularity) Let 1 < p < oo. If F £ L P (Q) and uq £ VF 2 ^ 1 with 8 u uq = on 

d£l, then there exists a unique solution of 
du 



dt 



(A + 7) it + F for a.e. t G (0, T), u(0) = u 



that satisfies 



dit 



dt 



+ II(^7MI^II»II?<c(||f|| ?+ kii; 2( ,_^J, 



where C is a positive constant independent of T and -F. 



(2.8) 



Inequality f|2.8|) was first established as Theorem 9.1 of [30] in Chapter IV, but the independency 
of C with respect to T is given later as Theorem 1.1 of [31] (see also Theorem 2.3 of |21j). 

Now we consider the well-posedness of the following linear parabolic system which contains 
(PUD as its special case. 

( 

d t y = Ay-V ■ (By) - V • (aSJz) + F in Q, 
dtz = Az — jz + 5y in Q, 

(2.9) 

d v y = 0, d v z = on S, 

y(x, 0) = y (x) z(x, 0) = z (x) x £ $7. 

Proposition 2.1. Let a G L°°(Q) and B G L°°(Q) N with B ■ v = on S. 

(%) ijfyo) ^0 € L 2 (f2) and -F G L 2 (Q), then system ()2.9[) admits a unique solution (y, z) G V^Q) x 
1/ 1 (Q) satisfying 



\y\\vHQ) + W z " 2 



VHQ) 



< (I2/0II + koll + ll^lll) ; 



(2.10) 



(it) Let2<p<oo. If F £ L p (Q), y G L P (Q) and z G W 2{l ~^ ) ' p (n) with d u z = on dQ, then 



system f|2.9j) admits a unique solution (y,z) G L P (Q) x Wp'^Q) satisfying 



\y\\p + 11*11^8,1(3) 



II. IIP <e c "( l«hl£+IUhll p ... + HFH p 



^ W 2{1 p'' p (U) 



(2.11) 



(Hi) Let p > N + 2. If F £ L°°(Q), y £ L°°(Q) and z £ W 1,P (Q) iwtfi <9„z = on 90, «ften 
system (|2.9p admits a solution (y,z) £ L°°(Q) x L°°(Q) satisfying 



y\\oo + \\z\\oo <e K (|yo|oo + lko||wA.?>(fi) + Halloo) , 



(2.12) 
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where k is given by (jl.9p and C = C(0). 



Proof. The existence of solution with respect to i/Oj-^o an d F in different function spaces can be 
deduced similarly as in [30] for which we omit here. We only show the required estimates with 
respect to time T. Since the proof for (|2.10p is similar to (|2.11j) . we need only to show (12. lip . 
Multiply the first equation of (|2,9p by \y\ p ~ 2 y and integrate over Q, to get 
d 



\y\l + J | Vy| 2 \yr 2 dx < C (l + \\a\\l + \\B\t) \y\ p p + C \\a\t \Vz\ p p + C\Ff p , (2.13) 



dt 



and in the same way, to get from the second equation of (|2.9p that 



d i ip , 



dx + \z\ p p <C\y\ p p . 



(2.14) 



Differentiate \Vz\ p with respect to t and take the second equation of (|2.9p into account again to 
obtain 

j t \Vz\ p p + J \Vz\ p ~ 2 \Az\ 2 dx<C \Vz\ p p + C (\y% + \z\ p p ) . (2.15) 
The inequalities (|2.13p - (|2.15p together with Gronwall's inequality lead to 



\y(;t)\ p p + \z(.,t)\l + \Vz(;t)\ p p <e 



Ck 



Vol 



(2.16) 



for all t € [0, T\. On the other hand, by the maximal regularity (|2.8|) for the second equation of 
PH . it follows that 

H41 + iia^ii; + \\4l < c (ikoii;^.,^^ + \\v% + ikiij;) , 

which together with (I2.16j) yields (12.11j) . 

Now we turn to the L°°-estimate ()2.12p . We first assume that yo G C (J7) and F € C (Q) . Let A 
be defined by ([HE]), and let {e~ tA } t>0 and {e"* (A+7) } t>0 be the analytic C -semigroups generated 



by —A and — (A + 7) in L p (f2), 1 < p < 00, respectively. Then the solution (y, z) of system (|2.9 
can be represented as follows 



e~ tA yo 



+ [ e- {t - s)A [-V ■ {By) -V ■ (aVz) + F](-,s)ds, 
Jo 



< t ^ A+ ^y(;s)ds. 



(2.17) 
(2.18) 



Take the norm of C ($7) on both sides of (|2.17p to get 

||y(-,i)Hc7(n) < h~ tA yo\\c(Ti) + 1 \\e- {t - s)A V-(By + aVz)(;. 



+ 



c(n) 



ds 



c n) 



ds. 



(2.19) 



To estimate (|2.19p . we first observe that the operator —A generates a bounded analytic semigroup 
on C (f2)([6]). It follows from the maximum principle that 



follc(n) ^ II2/0I 



c(n) 



(2.20) 
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and 



-(*-«) 



(2.21) 



ds 



**<•,•) c(n) <l|F(-.»)ll c(n) 

for any < s < t. Since p > N + 2, we can take e and a such that 

p-N-2 N 1 1 

< e < and — < a < e. 

2p 2p 2 p 

Then, with the help of (|2.3p . (|2.6p and the Holder inequality, we have, for any t € [0, T], that 

f e- (t - s)A V-( J By + aVz))(, S ) ds 

JO C^S2) 

< f (^ + 7 ) Q e- ( '" s)A (^ + aVz)(-, S 

JO 

< C [ (t-s)- a -^\(By + aVz)(;s)\ ds 

Jo 

< C(l + ||a|| 00 + || J B|| 00 )(||y|| p +||Vz|| p )r 
This together with (|2.16p gives 

f Q \e~^ s)A V ■ {By + oVz)) (•, s)|| da < e c * (\y \ p + \\z \\ w i, HQ) + . (2.22) 

By (p39]) - ([232l) . we obtain 

\\y\L < e CK (HyolL + INU P (n) + Halloo) • ( 2 - 23 ) 
Next, take the norm of W l,p (Q) on both sides of (|2. 18|) to get 

|p("j*)llwi,p(Q) ^ 



+<5 



(*-«)(^+7)y(. )S ) 



for any i E [0, T]. To estimate ()2.24p . we first notice that 



iy 1 -p(n) 



< e 



CT 



W 1 'P(Q) 



(2.24) 



(2.25) 



which can be obtained by the same energy method used in proving (|2TT6D . Let ± < a < 1 - \. By 
(|2.2p . (|2.5p and the Holder inequality, we have that for any t 6 [0,T], 



-(t-s)(A+ 



7) y(-, 



ds < C 



(A + 7 ) Q e- ( *- s)(A+7) y(-,s) 



ds 



< C / (t-s)- Q |y(-,s)l ds 



< C||y|LT- 



-a+l- 



This together with ()2.16p gives 



^)(^) y (,s)|| wip(n) ds < (|y | p + \\z \\ w i, m + \\F\\ p ) 



(2.26) 
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Finally, by ([2^) - (p^6l) and the Sobolev embedding W l ^{9) ^ C{9) for p > N, we get 

ML < e CK (llyolL + INIwn^n) + Halloo) ■ (2.27) 

To complete the proof, let us consider the general case that yo € L°°(f2) and F £ L°°(Q). This 
can be done by smoothing the data and density argument. Precisely, let {yonj^Li C C(S1) and 
{ F n}n=i c be such that Von -> 2/0 hi L 2 (0), F n ->• F in L 2 (Q) and |yonL < I2/0L » Halloo ^ 

Halloo • For each n, let (y n , z n ) be a solution of (|2.9p corresponding to yon, zo, F n , which satishes the 
inequalities (|2.10p and (|2.12p with (y, z) replaced by (y n , z n ). Thus, by the uniformly boundedness, 
we can extract subsequences of (y n , z n ) such that it converges to (y, z), which is a weak solution of 
(12. 9p corresponding to yo,ZQ and F. Moreover, y,z satisfy the inequality (|2.12p . □ 



3 Proof of Theorem 11.11 

To prove Theorem ll.il we first establish a global Carleman inequality for the adjoint system (jl.5p . 

Let uj' CC uj, that is, u/ C u, be a nonempty open subset. Then, there is a function (3 € C 2 (Q) 
such that (3(x) > for all x G Q, and (3\gQ = 0, |V/3(x)| > for all x 6 Q, \ uj' (see [HI Lemma 
1.1]). For A > 0, set 

e A/3 g A/3 _ e 2A||/3|| c(TT) 

ip = UT-ry a= t(r-t) ' (3 - 1} 



and 



7 (A) = e 2A ^llo ( n). (3.2) 



Lemma 3.1. Let fi £ L 2 (Q), i = 0, 1, ... , iV. T/ien i/iere exists a constant Xq = \q(£1,lj') > 1, 
suc/i f/iai for all A > Ao and s > 7(A) (T + T 2 ), 

J j [(sv) 1+d \Vz\ 2 + (s<p) 3+d \z\ 2 ] e 2sa dxdt 

(N 
j j (sip) d e 2sa \f \ 2 dxdt + J2 J j {s^) 2+d e 2sa \fi\ 2 dxdt 

+ J J (s<p) 3+d e 2sa \z\ 2 dxd?J (3.3) 
for all solutions z to the equation 

d t z-Az = f + j:li§- mQ, 
d v z = on E, 

z(x,0) = zq(x) x 6 fi, 

with zq £ L 2 (fi), where C = C(fi,a/), and 7(A) given by (|3.2p . 
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Essentially speaking, Lemma 13.11 has been proven in |23| (see also [22] ) but the explicit inde- 
pendency of the constant C with respect to T is shown in a similar way as in |16] and [18] , For 
notational simplicity in the sequel, we introduce 



and 



h(s, A; 



h(s, A; i 



(^) 3 |V^| 2 + M 5 |0| 2 



e 2sa dxdt, 



3 \a\2 



e 2sa dxdt. 



(3.4) 



(3.5) 



Lemma 3.2. There exists a positive constant \\ = C(Q, cj, cj')(1+ ||«|| 2 X) + H-BH^) satisfying 7(Ai) > 
Ai > 1 such that for any A > X\,s > 7(A) (T + T 2 ) and <fi T ,9 T € L 2 (Q,), the associated solution 
(cj), 9) to (jl.5p satisfies 



h(s, A; <f>) + I 2 (s, A; 0) < C\ II \ 8 (s<pfe 2sa \(f>\ 2 dxdt, 

J JQu, 



(3.6) 



where C\ = C\(fl, cj', cj). 



Proof. Applying Lemma 13.11 to the first equation of (jl.5p with d = 2 and the second one with 
d = 0, respectively, we obtain that there exist positive constants cq(Q,u') and \® satisfying 



7(A?) > A? = co(n,o/) (i + IH 
such that for all A > A? and s > 7 (A) (T + T 2 ), 



2 

00 



Wo >1 



/i(s,A;0) + / 2 (s,A;0) < ci 



( S c^) 5 |0| 2 + ( S <p) 



3 m|2 



e 2sa dxdt 



(3.7) 



(3.8) 



for all solutions (0, 9) to (j 1 . 5 [) with 4> T ,9 T G L 2 (f2), where and in what follows, the symbol Cf,i 
1,2,..., stand for some positive constants depending on Q,co' and cj. 

Next, let £ G C£°(ft) be such that £ = 1 in cj', £ = in O \ cj, < £ < 1 in cj, and 



A£ • r 1/2 G V£ • £" 1/2 € L°°(f2) Jv . 

The existence of such a function £ is easy to obtain (see, for instance [H]). Set 

n = (sipfe 2sa . 

Multiply the first equation of (11.51) by 6>??£ to get 

8 II (sipfe 2sa \9\ 2 £dxdt= II ri£9[-d t (f> - A<j> - BV<j)]dxdt 
JJq JJq 

= J J {7?£0 [-A0 + 7 + V • (a V<£)] +cj)9^{d t rf) + ri£9 (-A<f> - BV(j))} dxdt. 
Integration by parts gives 



(3.9) 



6 (sip) 3 e 2sa \9\ 2 idxdt = Ji, 



(3.10) 
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where 



Ji = [[ (j)9(dtV + 7r])Cdxdt, J 2 = [I 0V(??£) ■ VOdxdt, 

JJq JJq 

J 3 = - / / a0V(??£) • V(f>dxdt, Ja= 9V(i]^) ■ ^4>dxdt, 

J Jo JJq 



J s = - / / 9ri£BV(f)dxdt, J 6 = / / r/£V(9 • V(f)dxdt, 

JJq JJq 

J 7 = -JJ a\V4>\ 2 rj^dxdt. 
To estimate these integrals, we first observe by (|3.ip and (|3.9p that 

|^| < (s V fe 2sa ; \V( V 0\ < C(e /2 sV + i^)e 2sa . 
This together with Cauchy's inequality gives the estimation of Jj, i = 1, . . . , 6 as follows: 

Ji < eih(s, \ ; e) + -£- [f [M 3 + (sv) 7 } e 2sa |^| 2 £dxdt; (3.11) 



4ei JJq 

J2 <e 1 I 2 (s,X;0) + -^- JJ [(sip) 5 + X 2 (sip) 7 ] e 2sa \(b\ 2 Cdxdt; (3.12) 
Js < exh(s, A; 0) + ^ [(^) 3 + A 2 (^) 5 ] e 2sQ H 2 edxdt; (3.13) 

J 4 < e 1 J 2 (fl, A; 0) + II [(sip) 3 + X(sip) 5 ] e 2sa |V0| 2 (3.14) 
tel JJq 

h < eih(s, A; 9) + ° 3 ^ llo ° ^ ^ (sc/?) 3 e 2sa |V</>| 2 £eteeft; (3.15) 

Je < £ih(s, X ^) + ^jj (s^fe 2sa | V^| 2 £dxdt\ (3.16) 

«/ 7 < IklL f I (sip) 3 e 2sa \Vcf)\ 2 £dxdt, (3.17) 

where e\ is an arbitrary positive constant which will be determined later. 
The inequalities (l3TTil) - (l3T?] ) lead to 

J4 + J 5 + J 6 + J 7 < 3£iJ 2 ( S ,A;0) + -(l + ||a|| 2 o + ||B|| 2 X) ) 



X 2 (sip) 5 e 2sa \V<p\ 2 £dxdt. (3.18) 

Next, we estimate the integral on the right hand side of the inequality f|3. 18j) . Let 

rj = X 2 (s<pfe 2sa . 

Multiply the first equation of (jl.5p by fj£<j) and integrate over Q to obtain, by the integration by 
parts, that 

If X 2 (sipfe 2sa \V(j)\ 2 idxdt = S^J i , 

JJq i=8 
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where 



J s = -« // 4> 2 Zd t f]dxdt, J 9 = - II </>V(?)£) • V(j>dxdt, 
1 J Jq J Jo 

J w = II ( f)fj£B'V<f)dxdt, Jii = S / / 94>fj£dxdt. 



Since 

\d t fj\ < CX 2 (s^) 7 e 2sa , < C(e 1/2 A 2 sV 5 + A 3 sV 6 0e 2sQ , 

in the same way of estimating J1-J7, we can get for any £2 > that 

Js<C II X 2 (s<p) 7 e 2sa \<j)\ 2 tdxdt- (3.19) 



J9 < e 2 h{s,\-A) + ^-Jj Q [A 4 (^) 7 + A 6 (^) 9 ] e 2sa \<p\ 2 ^dxdt; (3.20) 
J10 < e 2 I 2 {s, A; 9) + ° JJ q X\s^) 7 e 2sa |0| 2 £dxdt; (3.21) 

Jli < e 2 h(s,X;(t>) + -^~ If A 4 (sy?) 7 e 2sQ |0| 2 ^dardt. (3.22) 



4e 2 77q 
Finally, we take 

8 5 5 

£1 = , and £9 = n o — x 

lOd' 10c lC5 (l + || a ||^ + 20c i 

to get, from (l3TT0]) - (^22l) . that 

A; 0) + J 2 (s, A; 9) < c 6 (l + ||a||^ + ||B||^) 2 1 1 X & (s^fe 2sa |0| 2 dxdt. 
Thus there is a positive constant 

7(Ai) > Ax = c 6 (l + \\a\t + \\B\t) > A? > 1 



such that for any A > Ai and s > 7(A) (T + T 2 ), the inequality (|3.6p holds, where A§ is given by 
(157711 ■ □ 



Proposition 3.1. There exist positive constants A and s such that, for all T > 0,(j) T ,9 T G L 2 (S1), 
i/te solution ((f>,9) of the system (jl.5p satisfies 

^(■,0)\ 2 2 + \9(-, 0)\ 2 2 <e CK II el™ \<P\ 2 dxdt, (3.23) 

where k is given by (|1.9p . 

Proof. By integration by parts, we observe that 

2 + |V^| 2 <(l + || J B|| 2 o )|0| 2 + 5 2 |0| 2 , (3.24) 



and 



^|^ + |V0| 2 + 2 7 |0| 2 <||a|| 2 oo |V^. (3.25) 
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Suppose first that {{aW^ > 1. Multiply (pT24"P by ||a||^ to get by (pT25j) that 

> 0. 



d 

Tit 



e C(l+||a||^+||B||^)t 



lallL \4>\l + \°\ 2 o 



Integrating above inequality over [0, t] for any t G (0, T] gives 



l|a|loolW-> u )l 2 
which implies that 



t)\l + \e(-,t)\l), 



M,o)\l + MM < ^[a+noii-.+iiBnDT+H.iL] + ]eM ^ 

for any t 6 (0, T]. The integration of (|3.26p on both sides over [T/4, 3T/4] leads to 



. )O )|2 + |0(.,o)|2<4 e c [( 1 +H 2 cc +ll B ll^+IHI 



Since 



it follows by (|3.6p that 



I 2 — rp 



(sip)~ 5 e~ 2sa , {s^)- 4 e- zsa < in O x 



\<p\ 2 + \6\ 2 dxdt. 



3-2sa 



T 3T 

7' T 



|0(-,O)|^ + |0(.,O)|^< 



where by taking A and s as 



2 C[(l+||q||g +||B||g )r+[|a[| no ]+^ 



A s (s^)V SQ |(/)| 2 dx^ 



Qu 



A = C l + ||a||^ + ||B||^ , S = C l + ||a||^ + ||B||^ (T + n, 



we get (13331) . 

Finally, if Halloo < 1, then 



d_ 



=C(l+||B||L)t 



H2 + H2 



> 0. 



(3.26) 



is a direct consequence of (|3.24|) and (|3.25|) . Thus, (|3.26|) verifies. In a similar argument as in the 
proof of || 

a||oo ^ 1; one can easily get (j3.23p . This completes the proof. n 
Proof of Theorem 11.11 Let s and A be such that the observability estimate (|3.23p and 

v (\) = e - A ^Hc(TT) < I (3.27) 

hold. Let e > and consider the following optimal control problem 

Minimize j J J |/| 2 e^dxdt + \ (\y(; T)| 2 + \z(-, T)\f) } 

subject to all / 6 L 2 (Q), where (y,z) is the solution of (|1.4p associated to /. The existence of an 
optimal pair (f e ,y £ ,z e ) to the above optimal control problem follows from the standard argument. 
By the Pontryagin maximum principle ([!]), 



(3.28) 
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(3.29) 



Here, ((j> e ,6 e ) is the solution of the adjoint system following: 

( -dtcPe = A0 e + BV4> £ + S9 e in Q, 

-d t 6 e = A8 E - 1 Q e - V ■ (aV&) in Q, 

d v (j) £ = Q,d v e e = Q onS, 

(0 £ ,8 £ )(x,T) = -l(y £ ,z e )(x,T) ieO, 

where (y £ ,z £ ) is the solution of CO} with / = f £ . By (fl~4"|) . (13381) . (I3T29I) . and Proposition [3J] it 
follows that 



12 gs-^dt + 1 ( | y (. >T )|* + |*(.,T)|* ) < ( |y |^ + \z \ z 2 



1 



(3.30) 



We can simply get from ()3.28|) and (|3.30p that the control function f £ satisfies 



„ 2 < e CK 

e 112 — e 



I2 + 1 -^o 1 i 



Next we show that f £ can be taken in L°°(Q). To this end, let r be a sufficiently small 



positive constant and let {tj}.^ 1 be a finite increasing sequence such that < tj < r,j = 
0,1,... , M, tm+1 = t. Let {pi}fl De another finite increasing sequence such that po = %,PM > 
(N + 2)/2 and, 

/AT \ / 1 1\ 1 

(3.31) 



N 



+ 1 



1 1 \ 1 

< -^ = 0,1 

K p i+1 ) 4 



? J * * * J 



M - 1. 



Set 



By (I3I|), 



1 - e 



2«ll/3|| 



C(S2) 



ao = T Q= t(T-t) 



«o < a < 



Of) 



< 0, 



1 + 77(A) 
where 77(A) is defined by (|3.27p . 

For each i, i = 0, 1, . . . , M, M + 1, define 

CiOM) = e (s+r * )a > £ (*,T-i), 
= e^) ao ^(x,T-t), 
G l (x,t) = [d t (e^ ao )] <j> £ (x,T-t), 

Hi(x,t) = [d t (e^ ao )]e £ (x,T-t), 

and 

a(x, t) = a(x, T — t),B(x,t) = B(x,T — t). 
Then for each i, (Q, q{) solves the following system: 

dtd - A& = BVd + Sgi + Gi in Q, 

d t Qi - Agi = -jQi - V • (oVCi) + #i in <5, 

d u Q = 0,d v Qi = on £, 

Ci(s,0) = 0, ft(s,0) =0 xe!l. 



(3.32) 
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Now we apply the LP-L q estimate to the above system. By the semigroup theory, the solution 
(Ci) f?i)> * = 1) 2, . . . , M + 1, of (|3,32p can be represented as 



C<(-,t) = 
£>*(•,*) = 

Firstly, by to (I3T33I) . we have 



-(t-s)A 



-(t-s)-A 



(•,s)ds, 



[-7<?i - V • (5VCi) + Hi] (•, s)ds. 



(3.33) 
(3.34) 



ft N f 1 1 

Ki(-,*)U = G y m(t-a) 2 ^-i «; 



Pi-i 



which can be estimated by Young's convolution inequality (see, e.g. p~3l p. 3]) as 

IKilU < G^II^UVCill^ + + 11^11^) 



m(t) 2 Vw-i Pi 







(3.35) 



where n = 1/[1 - (1/pi-i) + (l/p»)]- Similarly, applying (pl4"j) and 077]) with e = | to (pTM]) , we 
have 



\Qi(-,t)\ 



?t N I 1 1 ' 

G / m(f - s) 2 p - 
o 



-7ft + -H r t)(-,a)| pi _ 1 ds 



+G||a|| 00 / m(t-*) 2 ^ 3 |VCi(-,s)| n ^ 

JO 

which can also be estimated by Young's convolution inequality as 



ft|| w < Clllftll^ + 11^11^ 



T 



m(t) 2 p " 



+G||a|l||VO 



'Pi— i 



m(t) 



jV / 1 1 \ 3' 



(3.36) 



Owing to (|3.3ip . we also have 



m(t) 



df 



o 



< G 



(T + 1) n + T (^ +1 )(p« 



and 



AT / 1 1 \ 3 



m(t)V 2 Vfi-i n J A > l dt 



< G 



(T + 1) r i + T 



Secondly, we estimate the energy of solution (Q, qi) to get the following LP*- 1 -estimate 
This inequality together with (1335]) . ([336]) . (l337jl . and (1338]) gives 



w + II 8i \\ Pi < e ° K (ll^lU-! + ll-Hillpi-i 



(3.37) 



(3.38) 



(3.39) 



(3.40) 



1(3 



Since 

WGiWp^ < cruo-ill^ and imi^ < CTllft-ill^ 

it follows from dOU]) that 



IICillK + llftllp, <e c 'M HCi-iHp^ lift 



'Pi-i 



(3.41) 



(3.42) 



where Cj = Ci(Q, oj), i = 0, 1, . . . , M, are positive constants. The iteration inequality (]3.42p from 
to M implies that 

||Cm|| Pm + \\qm\\ pm < e CK CM Co list + IM| 2 ) • (3.43) 
By the definition of Co and qq, we obtain from (13.30p and (13.43P that 



IICmIL f + \\QM\L <e CK (||y || 2 + lkoll 2 ) 



(3.44) 



Finally, we apply LP M -maximal regularity for the first equation of (I3.32p for Cm+i to get 

\\d t (M+i\\ PM + ||ACa/+i|| Pa/ + ||Cm+ 
< c (ll^lL ||VCm+i|| Pm + \\qm+i\\ Pm + IIGm+iIU 

This, by taking into account of (13.4ip and (I3.39p . leads to 

Hm+i\\ w »a {q) < e CK (\\Cm\\ Pm + \\qm\\ Pm ) ■ 
Hence, by the imbedding inequality ([301 Lemma 3.3,Ch.II]) 

IICm+i|| C (q) < e c{1+T+ ^ \\Cm+i\\ w 2,i (q) , 
for p M > (N + 2)/2, and by (l3T4"4"D . we get 

||CM + l|loo<e CK (||yo|| 2 + Po|| 2 ). 



That is 



D (s+r)a 



<e^(|| yo || 2 + |No|| 2 ). 



which together with (I3.28P yields 

;-*(|-ij(A))+r(l+»»(A))]oy e 



< e 



Ck 



+ Ikolla), 



where 77(A) is given by (|3.27p . This gives, by choosing r small enough such that 



~ s o - vW + r(l + 77(A)) <0, 



that 



£ 11 00 



< e 



Ck 



+ NII 
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The above inequality enables us to extract a subsequences of f £ , still denoted by itself, such that 
f e — > / weakly in L 2 (Q), weakly* in L°°(Q) as e — s> 0. Let (y e ,z £ ) be the solution to the system 
associated to f £ . Then, by Proposition 12.11 we see that y £ and z £ are both bounded in V 1 (Q). 
Thus, there exist subsequences y £ and z £ , still denoted by themselves, such that 



V,Ze 



z weakly in V 1 (Q); strongly in L 2 (Q) 



for (y,z) G V 1 (Q) n C([0, T]; L 2 (Q,)), which is the weak solution of the system corresponding to 
/ G L°°(Q), and y(x, T) = and z(x, T) = almost everywhere in This completes the proof. □ 

4 Proof of Theorem 11.21 

Let (u,v) be a trajectory of the system (|1 .2j) with the initial value (Ho, Ho); which satisfies (|l,3p . 
Set it = II + y, v = v + z, yo = uq — Ho, zq = v o — Ho- Then, (y, z) solves the following parabolic 
system 

8 t y = Ay - X V • (yVv) - X V • ((u + y)Vz) + l u f in Q, 
= Az — ^fz + 5y in Q, 

(4.1) 

d u y = 0, (9^2 = on S, 

y(x, 0) = y (^) z{x, 0) = z (x) x G J). 

The local exact controllability of the system (jl.ip is equivalent to the local null controllability of 
the system (|4.ip . 

Let K = {rj € L°°(Q)\ < 1}. For each rj € K, we consider the following linearized system 



(4.2) 



8 t y = Ay - V • (By) - V • (a„Vz) + W in Q, 

dtz = Az — jz + 5y in Q, 

d u y = 0,d v z = on S, 

y(x,0)=y o (^) z(ai,0) =z (aO x G fi, 
where = x(u + 77) and i? = x^U. By fjl .3[) . we see that 

a v eL°°(Q), B G L°°(Q) N with B ■ j/ = on S. 

so system ()4.2p is casted into the exact framework of system (|1.4|) . Thus, we can apply Theorem 
11.11 to obtain that for each r) G K, there exists a pair ((y, z), f) which solves system (|4.2|) with 
y(x, T) = 0, z(x, T) = almost everywhere in fi. Here and in what follows, we denote by (y, z) the 
solution to system (|4.2p corresponding to / and r\ if there is no ambiguity. By (|1.8p . we see that 
the control functions are bounded as follows: 



<e CKO (|yo| 2 + ko| 2 ) 



(4.3) 
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where 



c [ 1 + T+- 



By (|2.12p of Proposition 12.11 and (|4.3p , we have the following estimate 



\y\\vHQ) + \\4vhq) + \\y\ 



+ z 



< e 



Ck 



(l2/oloo + ll- Z o|lv^ 1 '«JV(Q) 



(4.4) 



(4.5) 



For rj € K, define a multi-valued mapping A : K -»• 2 L 'W) by 



y e l 2 (Q) 



3/ satisfying (|4.3p such that (y, z) is 
the solution to (j4.2|) corresponding to 77 and /, 
and y(x, T) = z(x, T) = a.e. in 



We apply Kakutani's fixed-point theorem ([8j p. 7]) to the map A to prove Theorem 11.21 First, it 
is clear that K is a convex subset of L 2 (Q). By the argument above, we see that A(r/) is nonempty 
and convex for each rj E K. Moreover, by (|4.5|) . A(7?) is bounded in V 1 (Q) for each rj € K and 
hence A(rj) is a compact subset of L 2 (Q) by the Aubin-Lions lemma (j8j p. 17]). 

Next, we show that A is upper semi-continuous. To this purpose, let {??n}^Li be a sequence of 
functions in K such that r\ n —¥ rj strongly in L 2 (Q), and let y n £ A(r/ n ) for each n. Then, by the 
definition of A(rj n ), there exists f n for each n such that (y n , z n ) solves the following system 

d t y n = Ay n - V • (By n ) - V • (a^Vzn) + l u f n in Q, 

d t z n = Az n - 'yzn + 5y n in Q, 

d v y n = 0,d u z n = on E, 

y n (x, 0) = y (x) z n (x, 0) = z (x) x e J), 



(4.6) 



and y n (x,T) = z n (x,T) = for x € ^ almost everywhere. Moreover, the control f n satisfies 



fn I 



By (|4.7p and Proposition 12.11 we obtain 



< e 



Cko 



(I2/0I2 + ko| 2 ) • 



lynllv^Q) + ll Z nlly 2 (Q) — e " ( 1 2/0 1 2 + Ikolltf/i 



»(n) 



(4.7) 



(4.8) 



By (|4.7|) . (|4.8|) and applying the Aubin-Lions lemma again, we can get / S L°°(Q), y S V 1 (Q), 
z G V^ 2 (Q) and the subsequences of f n , y n , z n , still denoted by themselves, such that 

fn f weakly* in L°°(Q), and weakly in L 2 (Q); 

y n — ^ y weakly in V 1 (Q), and strongly in L 2 (Q); 

z n — > z weakly in V 2 {Q), and strongly in L 2 (0,T; H 1 ^)). 

Passing to the limit as n — > oo in (|4.6p . we get that (y, z) is a weak solution of (|4.6p corresponding 
to rj. We claim that that y G A(?]). Actually, let Y n = y n — y, Z n = z n — z, and F n = l UJ {f n — /). 
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Then (Y n ,Z n ) solves the following system 



d t Y n = AY n -V- (BY n ) 

~ v ' k?„VZ n + (a Vn - a v )Vz] + F n in Q, 

d t Z n = AZ n - -yZ n + 5Y n in Q, 

d v Y n = 0,d„Z n = onS, 

Y n (x, o) = o z n (x, o) = o x 6 n. 

Multiply the first equation of (|4,9p by Y n , and integrate over f2, to give 
r/ 



r// 



\Y n \l + \X7Y n \ 2 2 < CWBWllYnll + CWa.jllVZ^l 



In the same way to the second equation of ()4.9|) . we have 

d 



dt 



\Z n \i 2 + \VZ n \ z 2 + 1 \Z n \ A 2 <C\Y n \i 2 . 



Differentiate \V Z n \^ with respect to t to get, from the second equation of (|4.9p . that 

^ \VZ n \ 2 2 + \AZ n \ 2 2 + 7 \VZ n \l < C \Y n \ 2 2 . 
Since Halloo < C, it follows from (I4.10p -( l4.12p and Gronwall's lemma that 

\Y n (;t)\l + \Z n (;t)\l + \VZ n (;t)\l 

< e c(i+\\B\\i)T f r \ Vn _ v \2\ Vz f dx+ r Fn Y ndx ). 

\Jn Jo, J 

On the other hand, since (y, z) solves (|4,2p . by (ii) of Proposition 12. 1\ we get that 



which together with (j4.3|) implies 



(bo| p + l|^o|l w2(1 -i), P(n) + ||Wll p 



(O) 



(4.9) 



+C / \ Vn - if \Vz\ 2 dx + C [ F n Y n dx. (4.10) 



(4.11) 



(4.12) 



(4.13) 



(4.14) 



Since Wp A (Q) <->■ C^Q) for p > N + 2 ([U Lemma 3.3, Ch II]), it follows from (HTTP that 



|Vz|| c( Q )iV < c (j yo | p + IMI^a-i)^ 



(4.15) 



Since r/ n — > r/ strongly in L 2 (Q), Y n — > strongly in L 2 (Q), and F n — >■ weakly in L 2 (Q), thus, 
by (|4.15p . we see that the right hand side of (|4.13p tends to as n — > oo. Hence, |l^(-,i)| 2 — > 0, 
|Z n (-,i)| 2 — > for all t € [0,T]. Since y n (x,T) = z n (x,T) = in Q almost everywhere, we get that 
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y(x,T) = z(x,T) = in 17 almost everywhere, which implies that y G A(ry). This shows that A is 
upper semi-continuous. 

Now it remains to show that A(K) C K. By Proposition 12.11 for any y £ A(K), 

\\y\L<^ (boL + Nll^i-j>, (n) ) . 

where ci is a positive constant. Take 5 = e~ ClK ° such that if lunl + \\zn\\ , ri ii „ < 5 which is 

exactly ([l.lOp . then HyU^ < 1 and hence A(K) C iT. Therefore, the conditions of Kakutani's fixed 
point are satisfied, that is, if the initial data (uo,vo) satisfies (|1.10p . then there exists at least one 
fixed point y, which together with z, is the solution of (jl.ip corresponding with some control / and 
satisfies y(x, T) = and z(x, T) = for igl! almost everywhere. This completes the proof. □ 
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